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Анотація. Стаття присвячена розробці та аналізу апроксимаційно - ітераційних алгоритмів, заснованих 
на методі сіток та методі прямих, для розв'язання задачі оптимального керування еліптичною системою зі 
степеневою нелінійністю. Для числового розв'язання основної та спряженої крайових задач використано 
різницеві схеми другого порядку точності із застосуванням неявного методу простої ітерації. Обчислювальні 
схеми методу прямих для розв'язання зазначених еліптичних крайових задач реалізовані у поєднанні з методом 
стрільби для наближеного розв'язання крайових задач для відповідних систем звичайних диференціальних 
рівнянь, виникаючих у розглядуваній області після решітчастої апроксимації. Для мінімізації цільового 
функціонала використано відомі методи умовної оптимізації градієнтного типу (методи проєкції градієнта та 
умовного градієнта). Суть пропонованого апроксимаційно - ітераційного підходу полягає в заміні вихідної 
екстремальної задачі послідовністю сіткових задач, апроксимуючих її на сукупності сіток, що подрібнюються, 1 
застосуванні того чи іншого ітераційного методу градієнтного типу до кожної з «наближених» екстремальних 
задач. При цьому пропонується будувати лише декілька наближень до розв'язку кожної із «наближених» задач і 
приймати останнє з цих наближень, використовуючи кусково лінійну інтерполяцію за початкове наближення в 
ітераційному процесі для наступної «наближеної» задачі. Послідовність відповідних кусково лінійних 
інтерполянтів розглядається як послідовність наближень до розв'язку вихідної екстремальної задачі. У роботі 
обговорюються теоретичні основи такого комбінованого підходу, а також його переваги перед традиційними 
методами на прикладі розв'язання модельної задачі оптимального керування. 
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Вступ 

У даний час математичне 
моделювання є невід'ємною частиною 
скільки | - | небудь з великих  науково- 
технічних розробок. Математичні моделі 
оптимізації для систем із розподіленими 
параметрами -- це найбільш складний клас 
оптимізаційних задач, особливо для систем 
керування нелінійного | типу, значна 
прикладна важливість яких виявляється, 
наприклад, під час оптимізації процесів 
теплофізики, дифузії, фільтрації, теорії 
пружності, а також під час розв'язання 
обернених задач для рівнянь математичної 
фізики, що розглядаються у варіаційній 
постановці |1). Потреби математичного 


моделювання нелінійних оптимальних 
процесів зумовлюють необхідність 
розробки | та | розвитку | ефективних 


числових методів оптимального керування 
з використанням комп'ютерних систем. 
Слід відмітити, що загальної 
методики для розв'язання нелінійних 
оптимізаційних задач на сьогодні не існує 
через їх велику різноманітність залежно 
від виду правих частин диференціальних 
рівнянь різних об'єктів: стаціонарних, 
нестаціонарних, із запізненням, нелінійних 
за координатами та керуваннями, їх 
структури та критерію оптимальності. 
Через ці (фактори о для знаходження 
оптимального керування зазвичай 
використовують таку методику, що 
дозволяє визначати і прогнозувати загальні 
й | окремі результати | залежно від 
особливостей | об'єкта | керування та 
критерія оптимальності. 
Аналіз останніх досліджень |і 
публікацій 
Для систем керування еліптичного 
типу конструктивна теорія і методи 
отримання оптимальних розв'язків 
представлені, наприклад, у монографіях |2- 
6). Результати сучасних досліджень як 
теоретичного, так 1 прикладного характеру 
стосовно задач оптимального керування 
системами зі степеневими нелінійностями, 
містяться в роботах |7-10| та ін. Зокрема, у 


роботі |7| досліджено керовані нелінійні 
системи, що використовуються як 
наближення до моделі 0 Біна теорії 
надпровідності П-го типу у просторовому 
випадку, з якими тісно пов'язана система 
рівнянь пористого середовища. 
Встановлено властивість фінітності носія 
розв'язання задачі Коші для нелінійних 
нестаціонарних систем у тривимірному 
випадку. У |З| досліджено граничне 
оптимальне керування в контактній моделі 
зі степеневим тертям у антиплощини; 
отримано умови оптимальності залежно 
від степеневого показника і запропоновано 
обчислювальну схему, засновану на 
лінеаризації та методі нерухомої точки. У 
І9| розглянуто задачу оптимального 
керування степеневою  нелінійністю для 
диференціального рівняння Шредінгера, 
що виникає в нелінійній оптиці. На 
підставі глобальної коректності такого 
рівняння показано існування оптимального 
керування, доведено диференційованість 
цільового функціонала за Фреше та надано 
умови оптимальності першого порядку. У 
роботі (|10| для керованої системи, що 
описується рівнянням еліптичного типу зі 
степеневою  нелінійністю, де залежність 
функції стану від керування не 
диференційована | за  Гато, доведено 
розширену диференційованість зазначеної 
залежності та отримано необхідні умови 
оптимальності, що дає можливість 
охопити нові класи критеріїв 
оптимальності, зокрема, вивчити задачу 
граничного спостереження. За малих 
значень швидкості зростання нелінійності 
та розмірності області в таких системах 
похідна Гато існує, і виведення умов 
оптимальності (як необхідних, так і 
достатніх) | не | викликає труднощів, 
внаслідок о чого з можна користуватися 
відомими методами оптимізації систем, що 


описуються нелінійними  еліптичними 
рівняннями. 

Постановка проблеми 

У даній | роботі | розглядається 


наступна задача оптимального керування 
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еліптичною системою зі  степеневою 
нелінійністю Пп 1: мінімізувати 
функціонал 


Диуа |(усд-сО0)4М0я 


зи (Од -уд) 40 а) 

за умов 
-АУ(Д У (б є (дян), хе О; (2) 
уУ(хХ)20, хеГ, (3) 


де О с 1 - відкрита обмежена область з 
достатньо гладкою межею Г; 


2(х) є 1, (0), у(х) є 1, (0), /(9 є І,(0) - 
відомі функції для х-(х,1,)Є0;ц»0 - 
відома стала; и є и(х) - керування, 
визначене на множині 

С «Ди(х) є І, (0): и(хує У м.в.), (4) 


Ус і(О) - опукла замкнена множина; 
у(х)є у(х и) - функція стану системи, 
що відповідає керуванню и є(/ ; 
2 8 
2о--б у з Оператор Лапласа 
дх дбх, 
Нелінійності задач для станів у 
математичних | моделях оптимального 


керування зазначеного вигляду можуть 
бути обумовлені цікавими для практики 
випадками: наявністю стоків субстанції 
(наприклад, дифузія речовини в активних 
середовищах із поглинанням речовини за 
нелінійним законом, в яких дифундуюча 
речовина вступає в хімічні реакції із 
середовищем, що супроводжуються 
нелінійним стоком субстанції); процесами 
кліматизації крізь межу або дифузії крізь 
мембрану, що обмежує область; 
біохімічними процесами; під час аналізу 
ліній передач з витоком в електротехніці 
та ін. (121. 

У роботі |11) встановлено факт 
існування та  єдиності в  банаховому 
просторі У - Н(О)П1,(О0) узагальненого 
розв'язку у(х) є у(хуи) крайової задачі (2), 
(3) за будь-якого допустимого керування 
иє(/. З результатів |13| до того ж 
випливає, що функціонал (1) за умов (2)- 
(4) є сильно опуклим і двічі неперервно 
диференційованим за Гато на множині (/, 
причому його похідна Гато в кожній точці 
изи(х)є У може бути подана у вигляді 


Г(п)ефО)ч 2и(и(х) -у()), хєО, (5) 
де функція м(х)ем(хри) є узагальненим 
розв'язком спряженої системи 

-ДОх) ЗУ" ОМОЮ) а 2009-х), (6) 

хе(х,х)єЄ0; 
фмО)є-0, хегГ. (7) 
Таким чином, обчислення градієнта 
(и) функціонала (1) за умов (2) - (4) в 
точці и зи(х)є / потребує послідовного 
розв'язання відповідних еліптичних 
крайових задач з однорідними умовами 
Діріхле: нелінійної крайової задачі (2), (3) 
відносно функції стану системи 
УС) є у(ди) і спряженої лінійної крайової 
задачі (6), (7) відносно м(х) ем (хуи). 
Отже, сильно опуклий неперервний 
функціонал (и) за умов (2) - (4) на 
замкненій | опуклій | множині | (/ | з 
гільбертового простору | 1,(02) досягає 


нижньої грані /" -іпі Ли) у єдиній точці, 
тобто для будь-якого  с(х)є І,(0) та 
уУ()є (0) існує 
(по), Ух) 
задачі керування (1) - (4), такий, що 
Ли") з", де и'єй'(х)є С - оптимальне 


майже для всіх 


оптимальний | розв'язок 


керування, | У (х)-у(пи'), хеб - 


відповідна оптимальна траєкторія системи. 


Мета дослідження 

Зазначена робота присвячена 
розробці та аналізу нових апроксимаційно- 
ітераційних алгоритмів, заснованих на 
методі сіток та методі прямих, для 
розв'язання задачі оптимального керування 
еліптичною системою зі  степеневою 
нелінійністю (1) - (4). Для числового 
розв'язання еліптичних крайових задач (2), 
(3) та (6), (7) пропонується 
використовувати скінченно-різницеві 
схеми другого порядку точності із 
застосуванням неявного методу простої 
ітерації |П4) та обчислювальні схеми 
методу прямих |15)| у сполученні з методом 
стрільби для наближеного розв'язання 
крайових задач для систем звичайних 
диференціальних рівнянь | (ЗДР), що 
виникають в області О під час решітчастої 
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апроксимації. Для побудови  сіткових 
аналогів функціонала (1) використовується 
метод скінченних сум (квадратур), а для 
пошуку нижньої грані апроксимованого 
функціонала - відомі методи скінченно- 
вимірної умовної мінімізації. 


Виклад основного матеріалу 

Розглянемо питання про 
застосування проєкційно - ітераційного 
підходу до розв'язання задачі (1) - (4) 
оптимального | керування  еліптичною 
системою зі степеневою нелінійністю. При 
цьому в ролі проєкційних 
(апроксимаційних) методів пропонується 
використовувати методи сіток та прямих, а 
в ролі ітераційних методів мінімізації 
цільового функціонала - метод проєкції 
градієнта та метод умовного градієнта. 
Запропонований підхід укладається в 
загальну схему проєкційно - ітераційних 
методів розв'язання задач мінімізації з 
обмеженнями в Гільбертових просторах, 


теоретично обгрунтованих в |16|, і 
супроводжується детально описаними 
обчислювальними технологіями, що 
дозволяє довести ідею  проєкційно- 
ітераційного | підходу до успішного 
розрахунку. 

Проєкційний метод у гільбертовому 


просторі 
Нехай на деякій множині Ї/ дійсного 
сепарабельного гільбертового простору Н 
заданий обмежений знизу функціонал 
ЕК (и): 


ії Е(и)а Е'»-оо, (8) 
иєО 
Для наближеного розв'язання задачі 
мінімізації Е(и) на (/: 

К(иу-щії, йиєеС сн, (9) 
апроксимуємо функціонал К(и) 
послідовністю більш простих 
«наближених» функціоналів Е Ки), 
по1,2,., заданих відповідно на деяких 
множинах ії. дійсних  гільбертових 


просторів й. 


п 


Будемо вважати, що 
простори Н , ізоморфні підпросторам Н, 
основного простору Н 
(нене .спреас єн) і Ф - 


п 


лінійні неперервно оборотні оператори, які 


здійснюють взаємно однозначне 


відображення Н, на 8; Ф/ - оператори 


оборотного відображення Й, на Н,, 
причому Ф, рівномірно по п обмежені, 
тобто (|| Ф. || «С", п-1,2,.... 
випадку множини 0, с Н, мають вигляд 
О з(й є Й: й еФи,иє7 «СОН,), 
пої; СО. (10) 


припускати 


У цьому 


Будемо надалі, що 
функціонали Е(й,), па 1,2,... пов'язані з 
вихідним функціоналом  К(и) умовою 
близькості 
|Е (й) -ЕСФ,й )|«В,, Уй, є 0,, (11) 
де В, 0 при п-»ос, а множини 0, 
п 1,2,.. пов'язані з / умовою 
Миє 0 ЗІ до й, Є: о Фй, зи. (12) 
З умов (11) 1 (8), зокрема, виходить 
співвідношення 
Е"-В, «ЕсФй)-В, ак), ПЗ) 
ми, єй, зо Ом 
випливає обмеженість | знизу 
наближених  функціоналів 
Е (7, ) на відповідній множині О, - 


іі Р(й)оБ»-с, пеі,2,... 


й, Є, 


Будемо розглядати для кожного 


звідки 
кожного З 


п 1, 2,... задачу мінімізації Е (й,) на І. - 
Е(й,)-зіпі, йє0ОсН, (4 
п.1,2,... 
Збіжність в  Н  проєкційного методу 
розв'язання | задачі | (9) встановлює 
Теорема 1 (16). 

Нехай (/  - обмежена замкнена 
опукла множина з Ггільбертова простору 
Н, О, - множина виду (10) з гільбертова 
простору Н, (п а1,2,..), пов'язана з (/ 
умовою (12). Нехай виконується (8), 
функціонали | К(и) 1 Е (,) опуклі і 
неперервні на (/ 1 О, відповідно 
ро ої Нр ли, В 
близькості (11). Тоді за кожного п- І, 2,... 


задовольняють умову 


функціонал Е (й) досягає на Ї/ своєї 


нижньої грані і для будь-яких й, є М, (М. 
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- множина точок мінімуму Е (й,) на О,) 
шт Е (п)-Е"'. При шцьому будь-яка 


о 


послідовність (Ф/й/)., є мінімізуючою 


паї 
для К(и) на (/ і кожна її слабка гранична 
точка є точкою мінімуму К(и) на (/,ау 
разі єдиності точки мінімуму и' є / до неї 
слабко збігається вся послідовність 
(Фу), 


ЕК(и) о сильно 


Якщо, крім того, функціонал 


паї! 


опуклий на  (/, то 
йт Ф'й зи". 
пою 

Апроксимаційна схема, основана на 
методі сіток. 

Метод сіток (скінченних різниць) є 
найбільш вживаним числовим методом під 
час розв'язання крайових задач для 
еліптичних | диференціальних рівнянь 
завдяки його універсальності та 
ефективності. Універсальність полягає у 
можливості застосування цього методу як 
для лінійних, так і для нелінійних 
крайових задач, для задач різної 
розмірності та для областей складної 
(неканонічної) форми. 

Будемо вважати для простоти, що в задачі 
оптимального керування (1) - (4) область 
О с 12 є прямокутною, тобто 


Оз0)Г а хе(о,х): 0«х, 41, ас 1,2), 
і нелінійність у рівнянні стану (1) є 
слабкою, тобто існує стала с»0, така, що 
Зу (Х) «с, Ухе О. Вважатимемо також, 


що (/ - обмежена замкнена опукла 
множина з І, (0) , яка має вигляд 

ОС «Ди(х) є І, (0): а хи(х) «4, м.в.), (15) 
де д, і а, - задані сталі, -0«4 «4, «00, 

Для наближеного розв'язання задачі 

керування (1) (3), (15) введемо в О 
прямокутну рівномірну сітку вузлів з 
кроками Й, та ЛП, за напрямами х, та х, 
відповідно так, що Й, - І, / М, (а - 1,2), де 
М,»0 та М, » 0 -- цілі числа: 

хво: 
М,» 2-12). 


сукупність всіх 


в, -б, У, - Її", 
моїй, із 

Тут рю 1); о, 
внутрішніх вузлів; У, 


зо (а 


- межа сітки 0),: 


) Ю (й ї. хх» М-1 
Вес НО о З 


пгоранаю кова 171. замінимо 
інтеграли в (1) скінченними сумами за 
кубатурною формулою прямокутників, а 
рівняння стану (2)  -  різницевими 
рівняннями другого порядку точності |141. 
У результаті прийдемо до наступної задачі 
мінімізації сіткового функціонала 
м-1Мь-1 
/(ї)з У, уко Р є; Ру А, 
йо іззі 
М -1 М5-1 


оз зе Уа) ПЛ (16) 


зі іззі 


за умов 
-Лу,ь, а з ТИ» ХЕ; (7) 
Уь 70, хеу,, (13) 
йе0 «| йаи, |єН,: 4 хи, х4,, 
а 771,2,...М. 1 а-1,2) (19) 


у гільбертовому просторі Н, сіткових 
функцій й -и, ,), визначених на 6, Зі 


скалярним добутком елементів 
м-1 Мч 


(п, у), З У 3 ма БУЛ» Ми, уЄН, 
йзі іззі 
і нормою || д ||,,, 77 «/(7, й), - Тут 
ца ЗНО), У ЗУ, 
сь З м Же, а у, Я -убх, а 1 5 


дано. -1,М, -1 (а -1,2); 


М - різницевий оператор Лапласа: 


1 
Лу, що р на - ЗУ, ч Ук ї 


за ої 7 ду зу У, РРО хе о, 


Різницеву схему (17), (13) для фіксованого 
йє С з (19) можна записати у вигляді 
нелінійного операторного рівняння 


АузТЇчій (20) 
у просторі Н,, 
де Аре-лу, 


й із 


р хеФ,. 

Нескладно бачити, що оператор А 
диференційований за Гато в Н,, причому 
його похідна Гато А"(р) для будь-якого 


реН,: Р, 70, хєу, визначається за 
формулою 
пу 2 
АСРУУ З ЗЛУ, а ЯР Уа» Х Є 
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уу є Ні Уа 


Крім того, А"(р) є самоспряженим і 


-0, хеу,: 


додатньо визначеним лінійним оператором 
в Н,, | енергетично | еквівалентним 


різницевому оператору -Л: 
-р, А є Ар) є -и, ЛА 


де й2і й, 2ічс/у, У, мінімальне 
власне значення оператора Л. 
Для розв'язання нелінійного 


операторного рівняння (20) скористаємося 


неявним методом простої ітерації (14) 
о (ка) -(К) 


но 1 ад Хобі: 
т 
де візьмемо В--Л, рив і 


довільне початкове наближення - є Н,- 


Вочевидь, для знаходження 5") при 


рен 
розв'язати лінійне операторне рівняння 
ВУ"? -ф, февУ -к470- 7-й), 
або, що те ж саме, різницеву задачу 
Діріхле для рівняння Пуассона: 
Лу: с а -(1-т) Лу) з 
її» Уйі 


могу Даті) хєб; (21) 


у -0, хеу, Ке0,1,.. 


Оператор Л у схемі (21) можна вним 
одним із прямих методів, наприклад 
методом розділення змінних (14). За 
критерій закінчення ітераційного процесу 
(21) неявного методу простої ітерації 
сні) | ре ПА є, де 


заданому потрібно 


можна обрати такий: || У 


є» 0 - задана точність вн 

Після | отримання для вихідної 
нелінійної крайової задачі (2), (3) Її 
й й М ее 
сіткового розв'язку У -іу,,Ь хєб,, що 
відповідає даному керуванню  йєй/, 
можна розв'язувати  спряжену лінійну 
крайову задачу (6), (7) відносно функції 
мО) є Оки), х є О., Так само, як і раніше, 
замінімо еліптичне рівняння (6) на сітці 
5, 20, ОУ, відповідними  різницевими 
рівняннями другого порядку точності, в 
результаті чого прийдемо до сіткової 
крайової задачі 

2 
-Лм, із ЗУ ці ХУ, а), (22) 


хеф,; 


Мі З 0, хеу,» (23) 
де Л - різницевий оператор Лапласа, 

М Я ЧИО хо), У, Я Уа 
ну. ро Мі -1 (а-1,). 
Задача (22), (23) є системою лінійних 
алгебраїчних рівнянь відносно наближених 
значень Х/ - пу. . | функції м(х) е шм(хри) у 


й 
внутрішніх вузлах сітки 0,. 


о, бо 
ху), 


зь З 209 


152 


Розв'язок 
такої задачі існує, єдиний і стійкий по 
правій частині |14). Різницеву схему (22), 
(23) можна записати у вигляді лінійного 
операторного рівняння 

Тфе2о-7) 
у просторі Н,, де 


Ту о Лі, із зас 6 Чі ьо ХЕ оз 
Нескладно бачити, що Т є самоспряженим 
і додатньо визначеним лінійним 
оператором в Н, з межами спектра 


п, 7Ув По. 7: 76, 
де у, і У» відповідно мінімальне 1 
максимальне власні значення оператора Л : 


24 пл... 8. З 
іп'-о аа 


Ці властивості дозволяють 
застосовувати до розв'язання різницевої 
схеми (22), (23) різні ітераційні методи, в 
тому числі метод простої ітерації |14| 


(ен) (ю 2 
чи зЇ2 А, і тлу, із ра НА - 


а п) хе; (4) 


чат, хеуд Кеб,., 


де т22/(п,  т,), 
початкове наближення. За критерій 
закінчення процесу ітерацій (24) можна 
сен фо 

Їн є» є20. 


г (0) - 
у" єН, - довільне 


обрати такий: | / 

Після отримання для  спряженої 
крайової задачі (6), (7) її сіткового 
розв'язку ф - (м, , Б х є 9,, що відповідає 
керуванню йє (с Н,, можна визначити 
для даного йє( значення сіткового 
аналогу похідної Гато (5) на б,: 


ЛО) а ми, з У» х Є, (25) 


іі 
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і побудувати розрахункові формули деяких 
методів скінченновимірної умовної 
мінімізації для задачі пошуку нижньої 
грані сіткового функціонала (16) за умов 
П7) - (19). 

Розглянемо | метод наближеного 
розв'язання задачі оптимального 
керування (1) (3), (15), оснований на 
застосуванні до сіткової задачі (16) - (19) 
ітераційного методу проєкції градієнта. 
Такий метод приводить до побудови в Н, 
послідовності наближень (й), сі за 
формулами 

я 2 р (дО З-аЮТ(О 
ЗБ -аб ДО), (26) 
Кк 20,1... йФєЙ0, 
де Р. - оператор проєктування простору 
Н, намножину 0 ; //(й")) - похідна Гато 
функціонала | (16) в точці 
є Й, яка обчислюється за формулою 


сіткового 
по 


(25); а »0  - ітераційний параметр 
(кроковий множник) (17). З урахуванням 
зазначеного, формули (26) реалізуються 
так: 


(Ю Юр'усЮ 
и.-а ло )«а, 


й із 


бі 
(ен) (о Ю рсЮ 
цо'в34, йо о ЛО" )»4,, хе 


ць 150 
) 


й, -- 
о оди), ко0,. 


Нескладно перевірити, що 
сформульовані в |17| умови теореми про 
збіжність ітераційного методу проєкції 
градієнта виявляються для сіткової задачі 
оптимального керування (16) - (19) 
виконаними. 

Для о розв'язання  сіткової задачі 
керування (16) - (19), що апроксимує 
вихідну задачу оптимального керування 
() - (3), (15) на сітці 9,, можна також 


Ж 


К Кк . 
а ЛО), інакше, 


послуговуватися методом умовного 
градієнта (17), який призводить до 
побудови послідовності наближень 


(Оу? с О за формулами 
лен - а 


й Поч ой), 427) 

К 0,1, йФєй ; 
де 0«ха" «1 ітераційний параметр; 
лО е 0 - допоміжне наближення, яке 


вибирається з умови 


(ЛО, Оу, стіп (ЛИ), б, 


і яке для множини (19) виписується у 
явному вигляді: 


тут (Ю 
сф. Їй; 4120, 
й тут хеф,; 
а, Л(07)«0, 


б «аб К 0,1,... 


Збіжність методу умовного градієнта 
(27), (28) для сіткової задачі оптимального 
керування (16) - (19) випливає із загальної 
теореми про збіжність цього методу в 
гільбертовому просторі (17|. За критерії 
зупинки ітераційних процесів мінімізації 
(26) та (27), (28) можна взяти один або 
декілька з відомих критеріїв: 


СА и РИБ 
| Лу й72 7 МО)| єв; ДІЛУ, се, 


де є»0 - задана точність обчислень. 
Нескладно бачити, що метод сіток 
розв'язання задачі (1) - (3), (15) є за суттю 
методом  проєкційного типу, описаним 
вище. Дійсно, якщо позначити через 
ол по1,2,.. послідовність з сіток у 


(28) 


прямокутнику 0, що подрібнюються: 
Я а У, Кот у є: 


са 9? "а НА о -і, 2 2 (29) 
Ра М, 1,2 020 (с-1,2), 
дела ОМ що зі ІМ 


напрямами | х, | та | х, 


п 


ей НН В зе 1 


- кроки за 
відповідно, 
п, з (п, У, то в ролі задачі (9) 
(К(и)-з тв, иєИснН) туг 


задача мінімізації за умов (2), (3), (15) 
функціонала (1) ( К(и) - "(и)), заданого на 


виступає 


множині (15)  гільбертова | простору 
Н « 1.0, Т|, ароль наближених задач (14) 
СЕ (й, ) -зїпі, й, є0, с- Н, ), 
відповідно задачі мінімізації за умов (17)- 
(19) сіткових  функціоналів виду (16) 


відіграють 


(Ей )- м (й,)), заданих на множинах 


виду (19) з  гільбертових просторів 


Н, У Н,, сіткових функцій й, и и, і і 


хеф визначених на сітках (29), з 


я з 
нульовими значеннями на межі Ум 


(па12,...). 
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Дотримуючись  |17|,  підпростори 


Нн,еН, ізоморфні до просторів Н, 
визначимо як  гільбертові | простори 
кусково-сталих функцій 

мо? -1 му? Ж 

ща (п) О 
и (з У х ис Хі 00» хєО, (30) 
де ізз 
(п) 

де Ха, (2, хє С -  характеристична 
функція | комірки | ОО? -(ха(х,х,)єЮ: 


ді» 

5 | - 5 
рай «ай «хба"?, да1,2) сіткової області 
с, (п21,2,..). Функція (30) вочевидь в 
усіх точках такої комірки приймає 
однакові значення, рівні ие Легко 
бачити, що || и, |, 7- 11, ЇЇ я, » ТОБТО простори 


Н, і Н, ізометричні (18). Оператор Ф,, 
що відображає Н, на Н,, визначимо як 
оператор, який кожній функції и (х)є Н, 


виду (30) ставить у взаємно однозначну 


відповідність сіткову функцію 
й, -(и,, 1ЄН, » м. зи, Ох з хеф,? 
Ф'' - оператор, що здійснює оборотне 
відображення. Ясно, що || Ф, || - | Ф/ || -1 
(п с1,2,...). 


Як зазначалось вище, функціонал 
Ли) 1 7, (й,) сильно опуклі і неперервні 


на | замкнених | обмежених опуклих 
множинах М с 1.(О0) виду (15)1/, с Н, 


виду (19) відповідно (п- 1,2,..). Умова 


(12) для цих множин з урахуванням 


визначення оператора Ф/! виконується в 


силу сепарабельності простору (0), а 


здійсненність з умови близькості 0 (11) 
негайно випливає з твердження леми 9.2.2 
з П7). Таким чином, всі умови теореми 1 
виявляються виконаними, що гарантує 
збіжність | проєкційного  (різницевого) 
методу для задачі (1) - (3), (15). 


Апроксимаційна схема, основана на 
методі прямих 
Метод прямих належить до групи 
методів наближеного розв'язання 
крайових задач, головна ідея яких 
пов'язана зі зниженням розмірності задачі 
за рахунок апроксимації крайової задачі 
для | диференціального | рівняння | у 


частинних похідних системою 
диференціальних рівнянь, але з меншою 
кількістю неперервних змінних |151. 
Будемо застосовувати до розв'язання 
еліптичних крайових задач (2), (3) та (6), 
(7) при заданому иєї/ з (15) метод 
прямих, для чого побудуємо решітчасту 
область 0, з межею Г,, розбивши 


прямокутник О на М, --1 смуги прямими 
х, з х? з кроком Л, - І, / М, (М, » 0): 
Ом а оо б) єО: О«х «і, 
зе оо 3 
Межа ги 
відрізків прямих та точок: 
Гах «і, х, -0; Ок х хі, х, с; 


б.32); Сх 9). 12 М -Щ, (09) 
3. 


замінимо похідну --5- В 


2 
рівнянні (2) відповідним  різницевим 
виразом другого порядку точності |і 
отримаємо систему нелінійних ЗДР 


ЗЕ У, (237 2У, (0) У, а 0) 
ди з . 
? аруцним хни (88) 
о«х «Її, і, 21,М,-1 
з крайовими умовами 
у, (0) 20, у, (1) -0, і, -1,М, -1, 
у(2)20, у, (2220, Ох х хі, 


складається з таких двох 


На решітці О,, 


34) 


де у (цу), 1 є 
Аналогічно можна побудувати на 
решітці р. -О, Гу, лінійну крайову 
задачу для системи ЗДР, що апроксимує з 
другим порядком точності  спряжену 
еліптичну крайову задачу (6), (7): 


«Р, оз Мел СОМ Ск | 
9 (2) 
ах; 


п; 
Чу; Сх ум, ж) є 20, (2) сю з З; (35) 
бек 5-1 З 
.(0)-0, у, (1)-0, і, - 1, У, -1, 
у, (0) У, 1) й Й 36) 
уз (5) 20, у, (щ) 20, бхх хі, 


де щ, (0) я Ох), (що) є О,, 
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Зазначимо, що вибір напрямку 
прямих під час решітчастої апроксимації 
області С й основної та спряженої 
диференціальних задач на ній в загальному 
випадку залежить від вигляду нелінійності 
керованої системи і в даному випадку не 
має принципового значення. Відмітимо 
також, що для квазілінійних еліптичних 
задач розв'язуваність, стійкість, збіжність 
та оцінки похибки решітчастих схем 
досліджені у (|15|. Для наближеного 
розв'язування крайових задач (33), (34) та 
(35), (36) при заданому иє (/ з (15) можна 
скористатись методом стрільби (19, 20) з 
наступним числовим розв'язанням 
допоміжних задач Коші на сітці точок 


о, го єГ0, 21: хо нів, і 2 0,М,), 
де п -1/М, М»д0, за допомогою будь- 


якого відомого методу, зокрема методу 
Рунге - Кутта. 

Метод прямих розв'язання задачі 
керування (1) - (3), (15) є за суттю методом 
проєкційного з типу. Можна показати, 
користуючись наведеними вище 
результатами для методу сіток, що умови 
теореми 1 про збіжність проєкційного 
методу для задачі мінімізації функціонала 
(1) за умов (2), (3), (15) в цьому випадку 
також виконуються. 


Проєкційно - ітераційні схеми розв язання 
задачі керування 

Розглянемо питання про застосування 
до розв'язання задачі | оптимального 
керування (1) - (3), (15) проєкційно- 
ітераційного методу, основаного на методі 
проєкції градієнта. | Зазначимо, що 
стосовно до задачі опуклої оптимізації (8), 
(9) в гільбертовому просторі Н цей метод 
полягає в тому, що ітераційний метод 
проєкції градієнта використовується для 
розв'язання кожної із наближених задач 


п14) (з опуклими 
неперервнодиференційованими 
функціоналами | К(й,) | на опуклих 


множинах (/ ) в гільбертових просторах 


у 


Н 


п? 


наближень й? є0,, ГА о У 


причому будується лише декілька 


останнє 


п? 
з яких залучається ї до | формування 
початкового наближення в ітераційному 


процесі для Рено наближеної задачі в 
просторі 


па 


йно -Р, «о з во в'йО) 37 
реддв -1; по1,2,... 


о р ) У йо 
(й Ин -Ф,Ф,й , Па ї, зи т є б ). 
Умови слабкої та сильної зно вН 
послідовності | наближень | (Ф''й"? у", 


проєкційно-ітераційного методу (37) до 
точки мінімуму функціонала К(и) на / 
теоретично обгрунтовані в роботі |21| для 
випадка, коли кроковий множник б// » 0 


визначається з умови 
го кр (Ю г ру КТ) сен) | йо 2 
Е (й, 7-Е (й, 72 5 й, Й и, ур 


де 6.» 0 - параметр алгоритму. 
Застосування до кожної із 

наближених задач (14) з опуклими 

неперервнодиференційованими на опуклих 


множинах 0, функціоналами Е (й, ) 


ітераційного методу умовного градієнта за 
означеним проєкційно-ітераційним 
принципом приводить до побудови в Н- 
послідовності наближень |Ф, за а |з ЩО 


визначається формулами 
лану ще р сю тОЮ Ю 
й зйо 3 (иг'-й), (38) 
Кк «0,1,....К,-1; п-1,2,... 
йо ку) . й» 
-Ф Ф ЮР єй, ), 


Мі пчі 
де 0460) «1 - кроковий множник, який 
можна, наприклад обирати з умови 
Е9У- ЕТ З ей б | ЕФ (ці9У|, 
0 «8, «1 - параметр алгоритму. Збіжність 


проєкційно-ітераційного методу (38) для 
розв'язання загальної екстремальної задачі 
з обмеженнями в гільбертовому просторі 
обгрунтована | в | роботі | (16). | За 
запропонованою методологією можна 
перевірити виконуваність умов збіжності 
проєкційно-ітераційних процесів (37), (38) 
стосовно до задачі оптимального 
керування (1) - (3), (15) зі степеневою 
нелінійністю. 

Отже, для наближеного розв'язання 
задачі (1) - (3), (15) за допомогою будь-якої 
із наведених проєкційно-ітераційних схем 
(37), (38) слід до кожної із сіткових задач 
мінімізації (16) - (19) (або, що те ж саме, 
наближених задач (14)), починаючи із 
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задачі невеликої розмірності, застосувати 
відповідний ітераційний метод 


градієнтного типу та зоае на й, 


Чи )єб,, 


п 


З яких із 


лише декілька наближень й" 
Ка і,2,..к,» 
використанням кусково-сталої інтерполяції 


(30) обрати за початкове наближення 
она - (и) є 


Й 


останнє 


пчі? 
м) м 
и, Рі - У Х и | й (ж, 
ля 20 
до розв'язку наступної сіткової задачі (16)- 
(19) на сітці Критерієм зупинки 


(і) 
505), хесф, 


пої 


Пи 

проєкційно-ітераційних процесів (37), (38) 

може служити один або декілька з 
наступних критеріїв: 

2 дю з 

Іа, нь Є 

-К 
|, КО а» Л 


де є»0 - задана точність обчислень. 


ЦРУДАК СТР РАБ 


Ки 
(й я «є, 


Пи пч1 


Числова реалізація алгоритмів 

Для аналізу ефективності 
запропонованих | проєкційно-ітераційних 
схем було виконано їх програмну 
реалізацію та тестування на прикладі 
розв'язання декількох модельних задач 
оптимального керування виду (1) - (3), (15), 
одна з яких мала такі вхідні дані: 


о хеСО ої бах. 518-120), рії, 
4-1, 4,-3, с0-х 090, у -х1-х), 
ЛО) -5іп 5 5іп(1-х). Під час реалізації 


проєкційно-ітераційних алгоритмів 
сукупність | вкладених вдвічі | сіток 
о, (п1,2,...т) визначалася порядком 


МО «МО 25 первинного розбиття і для 
отримання наближеного розв'язку задачі з 
точністю є-107" виявилася складеною з 
то сіток з порядком дискретизації 
останньої М" «М? -80. Кількість К, 
наближень на  п-му кроці алгоритму 


визначалась як найменше ціле Кк, що 
задовольняло нерівність 


ло сае) 
| "ри (й, ) и "ри (й ) | чо Є, 7 

в якій є, » 0 вибиралось згідно з порядком 

сіткової апроксимації вихідної задачі 

керування; критерієм закінчення 

проєкційно-ітераційного процесу служила 


умова є, 2є. Зокрема, на початковому 


наближенні й? (І, 5), балу З і» ЩО 
хр а 


задавалось на першому кроці а РН 
ітераційного алгоритму (38), основаному на 

методі умовного градієнта, наближений 

розв'язок вихідної задачі було отримано за 
т 

УЮ, з 23343 3-15 ітерацій; час рахунку 
па! 


склав 3,2 хв. для функціонала (1) було 


отримано наближене значення 
Ли)я0,3075. Під час розв'язання тієї ж 
задачі керування звичайним  сітковим 


методом з використанням методу умовного 
градієнта на одній останній сітці з 
М, - М, 80, точністю обчислень є - 10 


початковим наближенням, рівним /1,2 в 
кожному вузлі такої сітки, наближений 
розв'язок задачі було знайдено за 23 
ітерації, що зажадало 8,30 хв. часу; для 
функціонала (1) було отримано наближене 


значення Ли) я 0,3098. Графічні 
зображення отриманих наближених 
розв'язків модельної задачі керування 


подані на рис. І, 2. 


10 


Рис. 1. Графік залежності керування И(Х,, Х, ) 


Рис. 2. Графік залежності функції (оо) 


системи 


стану 
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Висновки 

У роботі розглянуто результати 
застосування проєкційно-ітераційного 
підходу до розв'язування задачі 
оптимального керування системою, що 
описується диференціальним рівнянням 
другого порядку у частинних похідних 
еліптичного типу зі степеневою 
нелінійністю. Розроблено  проєкційно- 
ітераційні алгоритми, основані на методах 
скінченних різниць і прямих як методах 
проєкційного типу та ітераційних методах 
проєкції градієнта й умовного градієнта. 
Досліджено | умови апроксимації та 
збіжності сіткових алгоритмів, розглянуто 
різні стратегії подрібнення сіток і 
формування початкових наближень в 
проєкційно-ітераційних алгоритмах. 
Розроблено комп'ютерні засоби, що 
реалізують різні обчислювальні схеми 
проєкційно-ітераційних алгоритмів, 
проведено аналіз їх ефективності на 
прикладі розв'язання модельних задач. 

За отриманими у роботі результатами 
можна зробити такі основні висновки: 

е- в | проєкційно-ітераційних  алго- 
ритмах, основаних на методах сіток 
і прямих, не потрібно знати 
заздалегідь кінцевий порядок 
дискретизації області, визначення 
якого є практично складним; 

е в  проєкційно-ітераційних  алго- 
ритмах полегшується вибір 
початкового наближення, оскільки 
воно має обиратися лише для 
першої сіткової задачі невисокого 
порядку; початкове наближення для 
кожної із наступних сіткових задач 
однозначно визначається схемою 
методу. Якщо при цьому врахувати, 
що під час розв'язання прикладних 
нелінійних задач вибір початкового 
наближення не завжди вдається 
зробити вдало, то стає зрозумілим, 
що застосування | проєкційно- 
ітераційного підходу виявляється 
більш доцільним, оскільки він 
дозволяє з вже | при невеликому 
порядку сіткової задачі суттєво 
наблизитися до шуканого розв'язку. 
Остання перевага | набуває | ще 
більшого значення, якщо з ростом 


порядку дискретизації диференці- 

альної задачі обумовленість 

системи  сіткових рівнянь має 
тенденцію погіршуватися; 

- в  проєкційно-ітераційних  алго- 
ритмах зменшуються  обчислю- 
вальні | витрати | на | побудову 
наближень, значна частина яких 
будується для  сіткових задач 
невисоких порядків; 

е складність проєкційно-ітераційних 
алгоритмів несуттєво зростає в 
порівнянні зі складністю 
відповідних алгоритмів проєкцій- 
ного | типу, оскільки | перші 
відрізняються від останніх лише 
додаванням одного зовнішнього 
циклу з модифікації розбиття 
сіткової області; 

- до недоліків проєкційно- 
ітераційних | алгоритмів | можна 
віднести потребу у дещо більшому 
обсязі пам'яті комп'ютера при 
зберіганні з масивів ї змінних на 
проміжних етапах. 

Для подальших досліджень у цьому 
напрямку передбачається | розглянути 
можливість застосування регуляризуючих 
обчислювальних схем - за / проєкційно- 
ітераційним принципом для розв'язання 
некоректних задач оптимального 
керування системами з розподіленими 
параметрами. 
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